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Abstract 



We study the critical behaviors of period doublings in N (N = 2,3, 4, . . .) 
coupled inverted pendulums by varying the driving amplitude A and the cou- 
pling strength c. It is found that the critical behaviors depend on the range 
of coupling interaction. In the extreme long-range case of global coupling, 
in which each inverted pendulum is coupled to all the other ones with equal 
strength, the zero-coupling critical point and an infinity of critical line seg- 
ments constitute the same critical set in the A — c plane, independently of 
N. However, for any other nonglobal-coupling cases of shorter-range cou- 
plings, the structure of the critical set becomes different from that for the 
global-coupling case, because of a significant change in the stability diagram 
of periodic orbits born via period doublings. The critical scaling behaviors on 
the critical set are also found to be the same as those for the abstract system 
of the coupled one-dimensional maps. 
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I. INTRODUCTION 



The nonlinear dynamics of coupled nonlinear oscillators has attracted considerable at- 
tention in recent years. Such coupled oscillators are used to model many physical, chemical, 
and biological systems such as coupled p-n junctions Josephson-junction arrays [@], the 
charge- density waves ||, chemical-reaction systems j4|, and biological-oscillation systems 
||. They exhibit diverse bifurcations, multistability, chaos, pattern formation, and so on. 

The coupled nonlinear oscillators studied here are coupled inverted pendulums, consisting 
of ./V identical inverted pendulum coupled through some interaction mechanism. We first 
consider a constituent element of the coupled dynamical system, i.e., a single parametrically 
forced pendulum with a vertically-oscillating suspension point. It can be described by a 
normalized equation of motion ||, 

x = f(x, x, t) = —2irPQx — 2n(Q 2 — A cos 2irt) sin 2tix, (1) 

where a; is a normalized angle with range x G [0, 1), /3 is a normalized damping parameter, 
O is the normalized natural frequency of the unforced pendulum, and A is the normal- 
ized driving amplitude of the vertical oscillation of the suspension point, respectively. This 
parametrically forced pendulum has an "inverted" stationary state, corresponding to the 
vertically- up configuration with x = |. It is well known that as the parameter A is in- 
creased, the inverted pendulum undergoes a cascade of "resurrections," i.e., it becomes 
stabilized after its instability, destabilize again and so forth ad infinitum 0-lO|. Recently, 



we have studied bifurcations and transitions to chaos associated with such resurrections 
of the inverted pendulum |[L1|| . For each case of the resurrections, the stabilized inverted 
state exhibits an infinite sequence of period-doubling bifurcations accumulating at a period- 
doubling transition point A*, beyond which chaos sets in. Consequently, an infinite series of 
period-doubling transitions to chaos occur successively with increasing A. This is in contrast 
to the one-dimensional (ID) map [12], where only single period-doubling transition to chaos 
takes place. However, the critical scaling behaviors at each ith period-doubling transition 
point A* (i = 1,2,3,.. .) are the same as those for the ID map. 

In this paper we study the critical behaviors of period doublings in the system of iV (N = 
2,3,4, . . .) symmetrically coupled inverted pendulums by varying the driving amplitude A 
and the strength c of coupling between the inverted pendulums, and also compare them with 
those for the abstract system of the coupled ID maps [|T^,0]. The "coupling effect" of the 
strength and range of coupling on the critical behaviors are particularly investigated. Both 
the structure of the critical set and the critical scaling behaviors for the coupled inverted 
pendulums are found to be the same as those for the coupled ID maps found by one of us 



(Kim) and Kook |14 



This paper is organized as follows. We first introduce N symmetrically coupled inverted 
pendulums in Sec. [II] and discuss their dynamical symmetries and couplings. Bifurcations 
associated with stability of periodic orbits and Lyapunov exponents in the coupled inverted 
pendulums are also discussed in Sec. |TJ. We then investigate the critical behaviors of 
period doublings in the coupled inverted pendulums in Sec. [TV]. As in the single inverted 
pendulum [JXTJ , the coupled inverted pendulums undergo multiple period-doubling transitions 
to chaos [e.g., see Figs. 0(a), [](a) and ^](b) for the "stability trees" associated with the 
first, second, and third period-doubling transitions to chaos, respectively]. For each period- 
doubling transition to chaos, the critical behaviors vary depending on whether or not the 
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coupling is global. In the extreme long-range case of global coupling, the zero-coupling 
critical point with c = and an infinity of critical line segments lying on the line A = A* 
constitute the same critical set in the A — c plane, irrespectively of N. However, for any 
other nonglobal-coupling cases of shorter-range couplings, a significant change occurs in the 
stability diagram of 2 n -periodic (n = 0, 1, 2, . . .) orbits in the A — c plane, and consequently 
the structure of the critical set becomes different from that for the global-coupling case. It 
is also found that the critical scaling behaviors on the critical set are the same as those for 
the abstract system of the coupled ID maps [Q. Finally, a summary is given in Sec. [V|. 

II. SYMMETRIES AND COUPLINGS IN THE COUPLED INVERTED 

PENDULUMS 

In this section we introduce iV symmetrically coupled inverted pendulums and then 
discuss their symmetries and couplings. 

Consider N symmetrically coupled inverted pendulums with a periodic boundary condi- 
tion, 

%m fip^mi ^rai ^) 9{^mi -^m+lj • • • > •^m—l)j lj 2, ... , N. (2) 

Here the periodic boundary condition imposes x m (t) = x m+ ]y(t) for all to, the function 
f(x, x, t) is given in Eq. (H]), and g(xx, . . . , x^) is a coupling function, obeying the condition 

g(x, . . . ,x) — for all x. (3) 

The second-order differential equations (0) are reduced to a set of first-order differential 
equations, 

•^m = Vmi 

Vm = f{x m ,y m ,t) +g{x m ,x m+ i, . . . ,X m _i), TO = 1, 2, . . . , JV. (4b) 

Consider an initial orbit point z(0) [= (^i(O), . . . , z N (0))}, where Zi = (x iy yi) (i = 1, . . . , N). 
Then its Poincare maps can be computed by sampling the orbit points z(m) at the discrete 
time t = to (to = 1, 2, 3, . . .). We will call the transformation z(to) — > z(to+ 1) the Poincare 
map and write z(to + 1) = P(z(to)). 

The 2iV-dimensional Poincare map P has a cyclic permutation symmetry such that 

a^Paiz) = P(z) for all z, (5) 

where a is a cyclic permutation of z such that a(zi, z%, . . . , zn) = (z 2) . . . , z^, z\) and a -1 
is its inverse. The set of all fixed points of a forms a two-dimensional (2D) synchronization 
plane, on which 

Xl = ■ ■ ■ = x N , ?/! = ••. = y N . (6) 

It follows from Eq. (^) that the cyclic permutation o commutes with the Poincare map P, 
i.e., oP = Pa. Consequently, the 2D synchronization plane becomes invariant under P, i.e., 
if a point z lies on the 2D synchronization plane, then its image -P(z) also lies on it. An 
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orbit is called a(n) (in-phase) synchronous orbit if it lies on the 2D invariant synchronization 
plane, i,e, it satisfies 



x x (t) 



x N (t) = x*(t), y x {t) 



y N (t)=y*(t). 



(7) 



Otherwise it is called an (out-of-phase) asynchronous orbit. Here we study only the syn- 
chronous orbits. They can be easily found from the uncoupled inverted pendulum ([I]), 
because the coupling function g satisfies the condition (|3|). Note also that for these syn- 
chronous orbits, the Poincare map P also has the inversion symmetry such that 



SPS(z) = P(z) for all z, 



(8) 



where S(z) = —z. If a synchronous orbit {z(t)} of P is invariant under S, it is called a 
symmetric orbit. Otherwise, it is called an asymmetric orbit and has its "conjugate" orbit 
S{z(t)}. 

We now discuss the couplings between the inverted pendulums. Consider an element, say 
the mth element, in the N coupled inverted pendulums. Then the (m ± <5)th elements are 
called the 5th neighbors of the mth element. Here we consider the case where the coupling 
extends to the Kth [1 < K < y(^y^) f° r even (odd) N) neighbor(s) with equal strength. 
Hereafter we will call the number K the range of the coupling interaction. 

A general form of coupling for odd iV (N > 3) is given by 



K 



g(xi, ...,x N ) 



2K + 



- J2 [u(x 1+ i) - u(xi)} 



l=-K 



1 



2K 
K = h 



K 



- U ( X l+l) - u (%i) 



l=-K 

N - 1 



(9) 



where c is a coupling parameter and u is a function of one variable. Note that the coupling 
extends to the i^th neighbors with equal coupling strength, and the function g satisfies the 
condition (|3]). The extreme long-range interaction for K = is called a global coupling, 
for which the coupling function g becomes 



N 



g(%i, . . . , x N ) = — [u(x m ) ~ u(xi)] 

iV m=l 



N 



— J2u(x m ) -u(xx) 



(10) 



This is a kind of mean-field coupling, in which each inverted pendulum is coupled to all 
the other ones with equal coupling strength. All the other couplings with K < ^1 (e.g., 
nearest-neighbor coupling with K = 1) will be referred to as nonglobal couplings. The 
K = 1 case for iV = 3 corresponds to both the global coupling and the nearest-neighbor 
coupling. 

We next consider the case of even iV (N > 2). The form of coupling of Eq. holds 
for the cases of nonglobal couplings with K = 1, . . . , (N > 4). The global coupling 
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for K — y {N > 2) also has the form of Eq. (|10|), but it cannot have the form of Eq. @, 
because there exists only one farthest neighbor for K — y, unlike the case of odd N. The 
K = 1 case for N = 2 also corresponds to the nearest-neighbor coupling as well as to the 
global coupling, like the N = 3 case. 



III. STABILITY, BIFURCATIONS, AND LYAPUNOV EXPONENTS OF 

SYNCHRONOUS ORBITS 

In this section we first discuss stability of synchronous periodic orbits in the Poincare 
map P of the coupled inverted pendulums, using the Floquet theory JL5|. Bifurcations 
associated with the stability and Lyapunov exponents are then discussed. 

The stability analysis of an orbit in many-coupled inverted pendulums can be conve- 
niently carried out by Fourier-transforming with respect to the discrete space {m} 
Consider an orbit {x m (t) ; m = 1, . . . , N} of the N coupled inverted pendulums (0). The 
discrete spatial Fourier transform of the orbit is: 

1 N 

m=l 

j = 0,l,...,iV-l. (11) 

The Fourier transform £j(t) satisfies £j(t) = (,N-j(t) (* denotes complex conjugate), and the 
wavelength of a mode with index j is y for j < y and for j > y. Here ^ corresponds 
to the synchronous (Fourier) mode of the orbit, while all the other £j's with nonzero indices 
j correspond to the asynchronous (Fourier) modes. 

To determine the stability of a synchronous g-periodic orbit [x\(t) = ■ ■ ■ = 2tv(£) = x*(t) 
for all t and x*(t) = x*(t + q)), we consider an infinitesimal perturbation {5x m (t)} to the 
synchronous orbit, i.e., x m (t) = x*(t) +5x m (t) for m — 1, . . . , N. Linearizing the iV-coupled 
inverted pendulums (Q) at the synchronous orbit, we obtain: 

_ df(x*,x-,t) df(x*,x*,t) 

OX OX 

N 

+J2Gi(x*) 8x l+m - x , (12) 



i=i 



where 



dg(xx, ...,x N ) 
Gi(x) = - 



dxi 



(13) 



Xl=---=Xff=X 



Hereafter the functions G/'s will be called "reduced" coupling functions of g(x 1: . . . ,xjy). 
Let 5£j(t) be the Fourier transform of 5x m (t), i.e., 

1 N 

iV m=l 

j=0,l,...,iV-l. (14) 
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Here <5£o is the synchronous- mode perturbation, and all the other ££/s with nonzero indices j 
are the asynchronous-mode perturbations. Then the Fourier transform of Eq. (|l^) becomes: 



8$ 



^ + [ dx* 

^ /N n v 3 = 



dx* 



N 



1=1 



0,1, 



N-l. 



(15) 



Note that all the modes 5£/s become decoupled for the synchronous orbit. 
The equation (15) can also be put into the following form: 



8£j 
Sf)j 



J At) 



St* 

Sr)j 



0,1, 



iV — 1, 



(16) 



where 







V 



df(x*,x*,t) 
dx* 



N 



+ Y,Gi(x*)e 
i=i 



2m(l-l)j/N 



df(x*,x*,t) 
dx* 



(17) 



Note that each Jj is a g-periodic matrix, i.e., J jit) = Jj(t+q). Using the Floquet theory WE 



we study the stability of the synchronous g-periodic orbit against the jth-mode perturbation 
as follows. Let $>j(t) = (ty (t),(j)j(t)) be a fundamental solution matrix with $j(0) = /. 
Here (fij (t) and (jrp (t) are two independent solutions expressed in column vector forms, and 
I is the 2x2 unit matrix. Then a general solution of the g-periodic system has the following 
form 



6 Vj (t) 



= 0,1 N 



j = u, i,...,iv -1, (18) 
Substitution of Eq. ( |T8"D into Eq. ( p!6|) leads to an initial- value problem to determine $■,(£), 

= Jj(t)Qj(t), ^(0) = /. (19) 

Each 2x2 matrix Mj [= <&j(q)], which is obtained through integration of Eq. fliTf ) over the 
period q, determines the stability of the q-periodic synchronous orbit against the jth-mode 
perturbation. 

The characteristic equation of each matrix Mj (j = 0, 1, . . . , N — 1) is 



A?-trM,A, + detM,- = 0, 



(20) 



where trMj and detMj denote the trace and determinant of Mj, respectively. The eigen- 
values, Xj t i and Aj j2 , of Mj are called the Floquet (stability) multipliers, which characterize 
the stability of the synchronous g-periodic orbit against the j'th-mode perturbation. Since 
the j = case corresponds to the synchronous mode, the first pair of Floquet multipliers 
(-Vi, -V2) is called the pair of synchronous Floquet multipliers. On the other hand, all the 
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other pairs of Floquet multipliers are called the pairs of asynchronous Floquet multipliers, 
because all the other cases of j 7^ correspond to asynchronous modes. 
As shown in [ Dj , det Mj is given by 



det Mj = Jl tr Jjd ' = e - 27Tl3nq . (21) 

Note that all the matrices Mj's have the same constant Jacobian determinant (less than 
unity). Accordingly, each pair of Floquet multipliers (Xj t i, Xjp) (j = 0,1, . . . , N — 1) lies ei- 
ther on the circle of radius e _7r/3n,? or on the real axis in the complex plane. The synchronous 
periodic orbit is stable against the jth-mode perturbation when the pair of Floquet multipli- 
ers Xj,2) lies inside the unit circle in the complex plane. We first note that the Floquet 
multipliers never cross the unit circle in the complex plane and hence Hopf bifurcations do 
not occur. Consequently, the synchronous periodic orbit can lose its stability against the 
jth mode perturbation when a Floquet multiplier Xj decreases (increases) through —1(1) on 
the real axis. 

A more convenient real quantity Rj, called the residue and defined by 

1 + detM 7 - - trM 7 - . . 

Rj = —. 3 - r-t, j = 0, 1, . . . , N - 1, 22 

J 2 1 + detM, V ' 



was introduced in Ref. [18j to characterize stability of periodic orbits in 2D dissipative maps 
with constant Jacobian determinants. Here the first one Ro is associated with the stability 
against the synchronous-mode perturbation, and hence it may be called the synchronous 
residue. On the other hand, all the other ones Rj (j ^ 0) are called the asynchronous 
residues, because they are associated with the stability against the asynchronous-mode per- 
turbations. 

A synchronous periodic orbit is stable against the jth-mode perturbation when < Rj < 
1 (i.e., the pair of Floquet multipliers (A 3 - ( i, A^) lies inside the unit circle in the complex 
plane). When Rj decreases through (i.e., a Floquet multiplier Xj increases through 1), 
the synchronous periodic orbit loses its stability via saddle-node or pitchfork bifurcation 
(PFB). On the other hand, when Rj increases through 1 (i.e., a Floquet multiplier Xj 
decreases through —1), it becomes unstable via period-doubling bifurcation (PDB). We also 
note that a(n) synchronous (asynchronous) bifurcation takes place for j = (j ^ 0). For 
each case of the synchronous (asynchronous) PFB and PDB, two type of supercritical and 
subcritical bifurcations occur. For the supercritical case of the synchronous (asynchronous) 
PFB and PDB, the synchronous periodic orbit loses its stability and gives rise to the birth 
of a pair of new stable synchronous (asynchronous) orbits with the same period and a 
new stable synchronous (asynchronous) period-doubled orbit, respectively. However, for the 
subcritical case of the synchronous (asynchronous) PFB and PDB, the synchronous periodic 
orbit becomes unstable by absorbing a pair of unstable synchronous (asynchronous) orbits 
with the same period and an unstable synchronous (asynchronous) period-doubled orbit, 
respectively. (For more details on bifurcations, refer to Ref. ||19|| .) 

It follows from the condition (H) that the reduced coupling functions of Eq. fllBp satisfy 



N 

Y,Gi(x)=0. (23) 
i=i 
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Hence the matrix (|17]) for j = becomes 



Ja{t) — ( df(x*,±*,t) df{x*,±*,t) ] • (24) 

V dx* dx* / 

This is just the linearized Jacobian matrix for the case of the uncoupled inverted pendulum 



II]. Hence the synchronous residue Rq becomes the same as the residue of the uncoupled 



inverted pendulum, i.e., it depends only on the amplitude A. While there is no coupling 
effect on R , the coupling affects all the other asynchronous residues Rj (j ^ 0). 
In case of the global coupling of Eq. fllCf ), the reduced coupling functions become: 

where G(x) = -^u'(x). Substituting GVs into the second term of the (2,1) entry of the 
matrix Jj(t) of Eq. (|T7|) , we have: 

Vf^frV 27 ^- 1 )^ - / ° for J = 0, (2(]] 
2><W e ~ i -cu'(x) for j ± 0. W 



i=i 



Hence all the asynchronous residues Rj (j ^ 0) become the same, i.e., R\ = ■ ■ ■ = Rn-i- 
Consequently, there exist only two independent residues Rq and R\, independently of N. 
We next consider the non-global coupling of the form (§) and define 

G{x) = 2f^p[ u \ x )i (27) 
where 1 < K < (^y^) for even (odd) N larger than 3. Then we have 

-2KG(x) for I = 1, 
. G{x) for2</<l + K or 
G ' (A] -\ forN + l-K<l<N, ^ 

otherwise. 

Substituting the reduced coupling functions into the matrix Jj(t), the second term of the 
(2, 1) entry of Jj(t) becomes: 

N 

^(oOe 2 ^- 1 ^ = -S N (K,j)cu'(x), (29) 



i=i 



where 

S N (K, 3 ) = u — = 1 - (2A - +1)sinf . (30) 

Hence, unlike the global-coupling case, all the asynchronous residues vary depending on the 
coupling range K as well as on the mode number j. Since Sn(K,j) = Sn(K, N — j), the 
residues satisfy 



S 



Rj = R N -j, j = 0,1,..., AT- 1. (31) 

Thus it is sufficient to consider only the case of < j < f (^f±) for even (odd) N. 
Comparing the expression in Eq. ( |29| ) with that in Eq. (|26| ) for j ^ 0, one can easily see 
that they are the same except for the factor Sn{K,J). Consequently, making a change of 
the coupling parameter c — > Sn ^ k jy the residue Rj for the non-global coupling case of range 
K becomes the same as that for the global-coupling case. 

When the synchronous residue Rq of a synchronous periodic orbit increases through 1, 
the synchronous periodic orbit loses its stability via synchronous supercritical PDB, giving 
rise to the birth of a new synchronous period-doubled orbit. Here we are interested in such 
synchronous supercritical PDB's. Thus, for each mode with nonzero index j we consider 
a region in the A — c plane, in which the synchronous periodic orbit is stable against the 
perturbations of both modes with indices and j. This stable region is bounded by four 
bifurcation curves determined by the equations Rq = 0, 1 and Rj = 0, 1, and it will be 
denoted by Un- 

For the case of global coupling, those stable regions coincide, irrespectively of N and 
j, because all the asynchronous residues Rj's (j ^ 0) are the same, independently of N. 
The stable region for this global-coupling case will be denoted by Ug- Note that Ug itself 
is just the stability region of the synchronous periodic orbit, irrespectively of N, because 
the synchronous periodic orbit is stable against the perturbations of all synchronous and 
asynchronous modes in the region Ug- Thus the stability diagram of synchronous orbits of 
period 2 n (n = 0, 1, 2, . . .) in the A — c plane becomes the same, independently of N. 

However, the stable region Un varies depending on the coupling range K and the mode 
number j for the nonglobal-coupling cases, i.e., Un = Un(K,j). To find the stability region 
of a synchronous periodic orbit in the N coupled inverted pendulums with a given K, one 
may start with the stability region Ug for the global-coupling case. Rescaling the coupling 
parameter c by a scaling factor Sn ^ k j^ for each nonzero j, the stable region Ug is transformed 
into a stable region Un{K,j). Then the stability region of the synchronous periodic orbit is 
given by the intersection of all such stable regions Un's- 

Finally, we briefly discuss Lyapunov exponents of a synchronous orbit in the Poincare 
map P, characterizing the mean exponential rate of divergence of nearby orbits ||20|| . As 



shown in Eq. (0), all the synchronous and asychronous modes of a perturbation to a syn- 
chronous orbit becomes decoupled. Hence, each matrix Mj [= 3^(1)] with q = 1 determines 
the pair of Lyapunov exponents (<Jj,i, <Jj^) (j = 0, 1, N—l), characterizing the average ex- 
ponential rates of divergence of the jth mode perturbation, where o^i > Oj^- Since each Mj 
has the same constant Jacobian determinant (i.e., detMj = e~ 2nl3n ), each pair of Lyapunov 
exponents satisfies <Jj t i+<Jj t 2 = —2irj3Q. Note also that the first pair of synchronous Lyapunov 
exponents (o"o,i, 00,2) is just the pair of the Lyapunov exponents of the uncoupled inverted 
pendulum [ 1 1] , and the coupling affects only all the other pairs of asynchronous Lyapunov 



exponents (<x,- 1, Oj^) (j 7^ 0). Furthermore, all the pairs of the asynchronous Lyapunov expo- 
nents for the global-coupling case become the same [i.e., (01,1, cxi i2 ) = • • • = (crv-1,1? &n-i,2)], 
as in the case of the asynchronous residues. 
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IV. CRITICAL SCALING BEHAVIORS OF PERIOD DOUBLINGS 



In this section, by varying the two parameters A and c, we study the critical scaling 
behaviors of synchronous PDB's in the N symmetrically coupled inverted pendulums for 
P = 0.2 and Q = 0.1. It is found that the critical behaviors depend on the coupling 
range. In the global- coupling case, in which each inverted pendulum is coupled to all the 
other ones with equal coupling strength, the zero-coupling critical point and an infinity of 
critical line segments constitute the same critical set, independently of N. However, for 
any other nonglobal-coupling cases, the structure of the critical set becomes different from 
that for the global-coupling case, because of a significant change in the stability diagram of 
the synchronous 2 n -periodic orbits (n = 0, 1,2, . . .). The critical scaling behaviors on the 
critical set are found to be the same as those for the abstract system of the coupled ID maps 



14| . We thus consider separately two kinds of couplings, the global- and nonglobal-coupling 
cases. 



A. Global Coupling 

We first study the N globally-coupled inverted pendulums with the coupling function 
of the form (TTO) . As shown in Sec. ffTIl a synchronous periodic orbit is stable when all its 
residues Rj (j = 0, 1, ...,N — 1) defined in Eq. (|22|) lie between and 1 (i.e., < Rj < 1). 
Here Rq is the synchronous residue determining the stability against the synchronous-mode 
perturbation, while all the other ones Rj (j ^ 0) are the asynchronous residues determining 
the stability against the asynchronous-mode perturbations. For the globally-coupled case, 
all the asynchronous residues become the same, independently of j, and hence only one 
independent asynchronous residue (e.g., Ri) exists. Accordingly, the stability region of a 
synchronous periodic orbit becomes bounded by four bifurcation lines determined by the 
equations Rq — 0, 1 and R\ = 0, 1. Here the Rq = and 1 (_Ri = and 1) lines correspond 
to the synchronous (asynchronous) PFB and PDB lines, respectively. In such a way, we 
obtain the stability diagram of the synchronous 2 n -periodic orbits (n = 0,1,2,.. .) in the 
A — c plane. Note also that the stability diagram becomes the same, independently of N, 
because all the asynchronous residues Rj (j ^ 0) for each synchronous orbit are also the 
same, irrespectively of N. Consequently, the structure of the critical set and the critical 
behaviors for the global-coupling case become the same, independently of N. 

As an example, we consider a linearly coupled case in which the coupling function (|10D 

is 

1 N 

g(x l ,...,x N ) = c[—Y J x m -xi\- (32) 

iV m=l 



As in the uncoupled inverted pendulum [11], the coupled inverted pendulums exhibit mul- 
tiple period-doubling transitions to chaos. Here we study the first three period-doubling 
transitions to chaos. For each period-doubling transition to chaos, the zero-coupling critical 
point and an infinity of critical line segments constitute the critical set in the A — c plane. 
Three kinds of critical behaviors associated with the scaling of the coupling parameter c are 
found on the critical set, while the critical scaling behavior of the amplitude A is always the 



10 



same as that of the uncoupled inverted pendulum. Note that the structure of the critical 
set and the critical behaviors for the coupled inverted pendulums are found to be the same 
as those for the coupled ID maps 



Figure |T](a) shows the stability diagram of the synchronous orbits with low period q = 
1,2. The stable region of a synchronous orbit is bounded by its PDB and PFB lines. 
The horizontal (non-horizontal) solid and dashed boundary lines correspond to synchronous 
(asynchronous) PDB and PFB lines, respectively. Each bifurcation may be supercritical or 
subcritical. 

We first consider the bifurcations associated with stability of the synchronous inverted 
stationary point, corresponding to the vertically-up configuration (i.e., X\{b) — ■ ■ ■ — 5Cjv(t) = 
x*(t) = \ and yi(t) — ■ ■ • — yN(t) = U*(t) = 0). The inverted state is a symmetric one with 
respect to the inversion symmetry S. Its stability region is denoted by the IS in Fig. |l|(a). For 
the unforced case of A = 0, the inverted state is obviously unstable. However, when crossing 
the horizontal dashed boundary line of the IS, its first resurrection occurs, i.e., it becomes 
stabilized with birth of a pair of unstable synchronous asymmetric orbits with period 1 via 
subcritical PFB. (For more details on the resurrection of the inverted state, refer to Ref. |[L1|| .) 
This stabilized inverted state destabilizes again through asynchronous PDB and PFB when 
the nonhorizontal solid and dashed boundary curves are crossed, respectively. However, 
it becomes unstable via synchronous supercritical PDB when crossing the horizontal solid 
boundary line, and gives rise to the birth of a new synchronous orbit of period 2. This new 
synchronous 2-periodic orbit also is a symmetric one with respect to the inversion symmetry 
S, as shown in Fig. 0(b) and its stable region is denoted by the SP2 in Fig. 0(a). This 
synchronous symmetric orbit of period 2 loses its stability through asynchronous PFB's 
when crossing the non-horizontal dashed boundary curves. However, it becomes unstable 
via synchronous supercritical PFB when the horizonatl dashed boundary line is crossed, and 
consequently a pair of new stable synchronous orbits with the same period 2 appears. Note 
that the new pair of synchronous orbits is a conjugate pair of asymmteric orbits with respect 
to the inversion symmetry S, which is shown in Fig. 0(c). That is, the inversion symmetry 
is broken due to the symmetry-breaking PFB. The stable region of the asymmetirc orbits of 
period 2 is denoted by the ASP2 in Fig. 0(a). Each synchronous asymmetric 2-periodic orbit 
becomes unstable via synchronous supercritical PDB when the horizontal solid boundary 
line is crossed, and gives rise to the birth of a new synchronous asymmetric 4-periodic orbit. 
Here we are interested in such synchronous supercritical PDB's. 

Figure ^ shows the stability diagram of synchronous asymmetric orbits born by syn- 
chronous supercritical PDB's. Each synchronous asymmetric orbit of level n (period 2 n , 
n = 1,2,3, .. .) loses its stability at the horizontal solid boundary line of its stable region 
via synchronous supercritical PDB, and gives rise to the birth of a synchronous asymmet- 
ric period-doubled orbit of level n + 1. Such an infinite sequence ends at a finite value of 
A\ = 0.575 154 ■ • -, which is just the first period-doubling transition point of the uncoupled 
inverted pendulum ||11|| . Consequently, a synchronous quasiperiodic orbit, whose maximum 
synchronous Lyapunov exponent is zero (i.e., o"o,i — 0), exists on the A = A\ line. 

We examine the treelike structure of the stability diagram in Fig. |], which consists of 
an infinite pile of [/-shape regions and rectangular-shape regions. Note that the treelike 



structure is asymptotically the same as that in the coupled ID maps [14]. The [/-shape 



branching is repeated at one side of each [/-shape region, including the c = line segment. 
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The branching side will be referred to as the zero c side. However, the other side of each 
[/-shape region grows like a chimney without any further branchings (as an example, see 
the branch in Fig. 0(b) starting from the right side of the [/-shape region of the ASP2). 



As in the coupled ID maps JL4]], this rule governs the asymptotic behavior of the treelike 
structure. 

A sequence of connected stability regions with increasing period is called a "period- 
doubling route" []14[. There are two kinds of period-doubling routes. The sequence of the 



[/-shape regions with the zero c sides converges to the zero-coupling point c = on the 
A = A\ line. It will be referred to as the U route. On the other hand, a sequence of 
rectangular regions in each chimney converges to a critical line segment on the A = A\ 
line. For example, the sequence of the rectangular regions in Fig. [|(b) converges to a 
critical line segment joining the left end point q (= 3.427 742 • • •) and the right end point 
c r (= 4.796 277- • •) on the A = A\ line. This kind of route will be called a C route. Note 
that there are infinitely many C routes, while the U route converging to the zero-coupling 
critical point (A*,0) is unique. Hence, an infinite number of critical line segments, together 
with the zero-coupling critical point, constitute the critical set. 

We now study the critical behaviors on the critical set. First, consider the case of the U 
route ending at the zero-coupling critical point. We follow the synchronous orbits of period 
q = 2 n up to level n = 8 in the U route, and obtain a self-similar sequence of parameters 
(A n , c n ), at which each orbit of level n has some given synchronous and asynchronous residues 
Ro and R± (= R2 = ■ ■ ■ = Rn-i) (e.g., Rq = 1 and -Ri = 0). Then the sequence {(A n , c n )} 
converges geometrically to the zero-coupling critical point (Al,0). As in the uncoupled 
inverted pendulum fll|| , the sequence {A n } obeys a scaling law, 



AA n ~ 5~ n for large n, (33) 
where AA n = A n — A n _i and 5 ~ 4.67. The value of the scaling factor 5 agrees well with 



the Feigenbaum constant (= 4.669 ■ • ■) of the ID map |12j| . We also note that the sequence 
{c n } obeys a scaling law, 

Ac n ~ fi~ n for large n, (34) 

where Ac n = c n — c n _i The sequence of the scaling factor {fi n } (= Ac n / Ac n+ i) of level n is 
listed in the second column of Table [I] and converges to a constant \i (~ —2.5), which agrees 
well with the coupling-parameter scaling factor a (= —2.502 • • •) of the coupled ID maps 
near the zero-coupling critical point |TJ]]. It has been also shown in [TJj] that the scaling 
factor a is just the largest relevant "coupling eigenvalue" (CE) of the zero-coupling fixed 
map of the renormalization transformation for the case of the coupled ID maps. 

We also study the coupling effect on the asynchronous residue i?i n of the synchronous 
orbit of period 2 n near the zero-coupling critical point (A*,0). Figure ^| shows three plots 
of R\ >n {A\, c) versus c for n = 5,6 and 7. For c = 0, Ri tTl converges to a constant R\ 
(= 1.300 59. . .), called the critical asynchronous residue, as n — > 00. However, when c is 

nonzero R^ n diverges as n — » 00, i.e., its slope S n (= ' ! "" 



) at the zero-coupling 

(A*,0) 



dc 

critical point diverges as n — > 00. 

As in the scaling for the coupling parameter c, the sequence {S n } also obeys a scaling 
law, 



12 



S n ~ v n for large n. (35) 

The scaling factor v n (= S n+ i/S n ) of level n is listed in the third column of Table p] and 
converges to a constant v (~ —2.5) as n —>■ oo. Note also that the value of v agrees well 
with that of the largest relevant CE a of the zero-coupling fixed map. 

We next consider the cases of C routes, each of which converges to a critical line segment. 
Two kinds of additional critical behaviors are found at each critical line segment; the one 
critical behavior exists at both ends and the other critical behavior exists at interior points. 
In each C route, there are two kinds of self-similar sequences of parameters (A n , c n ), at which 
each synchronous orbit of level n has some given synchronous and asynchronous residues Rq 
and R\] the one converges to the left end point of the critical line segment and the other 
converges to the right end point. As an example, consider the C route in Fig. 0(b), which 
converges to the critical line segment with two ends (A\,c{) and [A\,c r ). We follow, in the 
C route, two self-similar sequences of parameters, one converging to the left end and the 
other converging to the right end. In both cases, the sequence {A n } converges geometrically 
to its accumulation value A\ with the ID scaling factor 5 (~ 4.67) like the case of the U 
route. The sequences {c n } for both cases also obey the scaling law, 

Ac n ~ /i _n for large n, (36) 

where Ac n = c n — c n _i. The sequence of the scaling factor \i n (= Ac n /Ac„ + i) of level n is 
listed in Table [TI|, and converges to its limit value \x (—2). We also note that the value of 
\i agrees well with that of the coupling-parameter scaling factor [y = 2) of the coupled ID 
maps near both ends of each critical line segment [|14j]. It has been also shown in [JTJJ] that 



the scaling factor v (= 2) is just the only relevant CE of a nonzero-coupling fixed map of 
the renormalization transformation for the case of the coupled ID maps. 

Figure 0(a) shows the behavior of the asynchronous residue Ri <n (Al, c) of the synchronous 
orbit of period 2 n (n = 5,6,7) near the critical line segment in Fig. 0(b). Magnified views 
near the both ends q and c r are also given in Figs. 0(b) and 0(c), respectively. For c = q 
and c r , Ri tn converges to a critical asynchronous residue R{ (= 0) as n — > oo, which is 
different from that for the zero-coupling case. The sequence of the slope S n of Ri >n at both 
ends obeys well the scaling law, 

S n ~ v n for large n. (37) 

The two sequences of the scaling factors v n (= S n+ i/S n ) of level n at both ends are listed 
in Table [Til], and converge to their limit values v ~ 2, which agrees well with the only CE 
[y = 2) of the nonzero-coupling fixed map governing the critical behavior at both ends 
for the case of the coupled ID maps. However, for any fixed value of c inside the critical 
line segment, R^ n converges to a critical asynchronous residue R\ (= 0.5) as n — > oo [see 
Fig. 0(a)]. This case of R\ = 0.5 corresponds to the superstable case of = (A^: the 
critical asynchronous Floquet multiplier) for the coupled ID maps [Q, because Eq. ( |22"D of 
R for the case of 2D maps reduces to the equation of R = 0.5 x (1 — A) for the case of ID 
maps. We also note that as in the case of the coupled ID maps, there exists no scaling factor 
of the coupling parameter inside the critical line segemnt, and hence the coupling parameter 
becomes an irrelevant one at interior critical points. Thus, the critical behavior inside the 
critical line segment becomes the same as that of the uncoupled inverted pendulum (i.e., 
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that of the ID map), which will be discussed in more details below. This kind of ID-like 
critical behavior was found to be governed by another nonzero-coupling fixed map with no 
relevant CE for the case of the coupled ID maps (14[ . 

There exists a synchronous quasiperiodic orbit on the A = A\ line. As mentioned in 
Sec. PH , its synchronous Lyapunov exponents are the same as the Lyapunov exponents of the 
uncoupled inverted pendulum, i.e., <7o,i = and 00,2 = — 27T/3S1. The coupling affects only the 
pair of asynchronous Lyapunov exponents (o"i,i, 01,2) [= (°"2,i, 02,2) = ■ • • = (ojv-i.i, 0^-1,2)]- 
The maximum asynchronous Lyapunov exponent 0^ near the critical line segment in 
Fig. ^|(b) is shown in Fig. [5]. Inside the critical line segment (q < c < c r ), the syn- 
chronous quasiperiodic orbit on the synchronization plane becomes a synchronous attractor 
with 0i 5 i < 0. Since the dynamics on the synchronous attractor is the same as that of the 
uncoupled inverted pendulum, the critical maps at interior points exhibit essentially ID-like 
critical behaviors, because the critical behavior of the uncoupled inverted pendulum is the 
same as that of the ID maps ||11||. However, as the coupling parameter c passes through q 



and c r , the maximum asynchronous Lyapunov exponent 01.1 of the synchronous quasiperi- 
odic orbit increases from zero. Consequently, the synchronous quasiperiodic orbit ceases 
to be an attractor outside the critical line segment, and the system of the coupled inverted 
pendulums is asymptotically attracted to another synchronous rotational attractor of period 
1. 

What happens beyond the first period-doubling transition point A\ is also interesting. 
As in the uncoupled ID inverted pendulum ||11|| , with increasing the amplitute A further 



from A = A*, the unstable inverted state undergoes a cascade of resurrections, i.e., it will 
restabilize after it loses its stability, destabilize again, and so forth ad infinitum. For each case 
of the resurrections, an infinite sequence of PDB's leading to chaos follows. Consequently, 
the coupled inverted pendulums exhibit multiple period-doubling transitions to chaos. 

As the first example, we consider the second period-doubling transition to chaos. Figure 
|5](a) shows the second stability diagram of the synchronous inverted stationary point and 
asymmetric orbits of level n (period 2 n , n = 0,1,2,3) in the A — c plane. When cross- 
ing the horizontal solid boundary line of its stability region IS, the unstable inverted state 
restabilizes with birth of a new unstable synchronous symmetric orbit of period 2 via syn- 
chronous subcritical PDB. This is the second resurrection of the inverted state. However, 
when the horizontal dashed boundary line is crossed, the stabilized inverted state becomes 
unstable via synchronous supercritical PFB, which results in the birth of a conjugate pair 
of synchronous asymmetric orbits with period 1. Then each synchronous asymmetric orbit 
of level n becomes unstable at the horizontal solid boundary line of its stability region via 
synchronous supercritical PDB, and gives rise to the birth of a synchronous asymmetric 
period-doubled orbit of level n + 1. Such an infinite sequence terminates at a finite value 
of A^(= 3.829 784 • • •), which is the second period-doubling transition point of the uncou- 
pled inverted pendulum [|TTJ. Note that the treelike structure of the stability diagram in 



Fig. |6](a) is essentially the same as that in Fig. 0(a). Hence, the critical set also consists of 
the zero-coupling critical point and an infinite number of critical line segments, as in the 
first period-doubling transition case. In order to study the critical behaviors on the critical 
set, we follow the synchronous asymmetric orbits up to level n = 7 in the U route and in 
the rightmost C route. It is found that the critical behaviors are the same as those for the 
first period-doubling transition case. That is, there exist three kinds of critical behaviors at 
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the zero-coupling critical point, both ends of each critical line segment and interior points. 

As the second example, we also consider the third period-doubling transition to chaos. 
The third stability diagram of the synchronous orbits with q = 1, 2, 4, 8 is shown in Fig. |^(b). 
A synchronous subcritical PFB occurs when crossing the horizontal dashed boundary line 
of the IS. Consequently, the unstable inverted state restabilizes with birth of a pair of 
unstable orbits with period 1. This is the third resurrection of the inverted state. However, 
the stabilized inverted state becomes unstable via synchronous supercritical PDB when 
the horizontal solid boundary line of the IS is crossed, and gives rise to the birth of a 
symmetric 2-periodic orbit. The subsequent bifurcation behaviors are the same as those for 
the first period-doubling transition to chaos. That is, a third infinite sequence of synchronous 
supercritical PDB's follows and ends at a finite value ^3(= 10.675 090 • • •), which is the third 
period-doubling transition point of the uncoupled inverted pendulum [[11] . Note also that the 



treelike structure of the third stability diagram is essentially the same as that in Fig. 0(a). 
Hence, the critical set is composed of the zero-coupling critical point and an infinity of 
critical line segments. Furthermore, the critical behaviors on the critical set are found to be 
the same as those for the first period-doubling transition case. 



In addition to the linear-coupling case (|32|), we have also studied two other nonlinear- 
coupling cases, 

1 N 

g(x 1 ,...,x N )=c[—'£ l x? n -Xi], n = 2,3. (38) 

iV m=l 

First stability diagrams of the synchronous orbits for the cases of the quadratic and cubic 
couplings are shown in Fig. ^(a) and 0(b), respectively. Their treelike structures are essen- 
tially the same as that in Fig. 0(a). Hence, the zero-coupling critical point and an infinite 
number of critical line segments constitute the critical set for each nonlinear-coupling case. 
Moreover, the critical behaviors for these nonlinear-coupling cases are also found to be the 
same as those for the linear-coupling case. 



B. Nonglobal Coupling 

Here we study the nonglobal-coupling cases with the coupling range K < y(^r^) for 
even (odd) N. The structure of the critical set becomes different from that for the global- 
coupling case, because of a significant change in the stability diagram of the synchronous 
orbits with period 2 n (n = 0, 1, 2, . . .), as will be seen below. 

As an example, we consider a linearly-coupled, nearest-neighbor coupling case with K = 
1, in which the coupling function is 

g(x\, . . . , xn) = ~(x2 + xn — 2xi) for N > 3. (39) 
3 



As shown in Sec. |H|, the stable region Un, in which a synchronous orbit is stable against 
the perturbations of both modes with indices and j 0), varies depending on the mode 
number j, because the asynchronous residue Rj (j ^ 0) depends on j. To find the stability 
region of the synchronous orbit, one can start with the stability region Uq for the global- 
coupling case. Rescaling the coupling parameter c by a scaling factor 1/5*^(1, j) [Sn(K,j) 
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is given in Eq. (pOQI, the stable region Uq is transformed into a stable region f/jv ( 1 , J ) • Then 
the stability region of the synchronous orbit is given by the intersection of all such stable 
regions U^s. 

As an example, we consider the case with N = 4. Figure || shows the stability regions of 
the synchronous asymmetric 2 n -periodic (n = 1,2,3,4) orbits. Note that the scaling factor 

1 linn lfn "v n mi t-v-i i -i i-vi t rn 1 1 1 /~\ 



s . , ! ( .v has its minimum value | at j — 2. However, for each synchronous orbit, U^(l, 2) itself 
cannot be the stability region, because bifurcation curves of different modes with nonzero 
indices intersect one another. We now examine the structure of the stability diagram in 
Fig. starting from the left side of the stability region of the synchronous asymmetric orbit 
of level 1 (n = 1). For the case of level 2 (n = 2), the zero c side of f/ 4 (l,2) including 
a c = line segment remains unchanged, whereas the other side becomes flattened by the 
bifurcation curve of the asynchronous mode with j = 1 Due to the successive flattening 
with increasing level n, a significant change in the stability diagram occurs. Of the infinite 
number of period-doubling routes for the global-coupling case, only the U route ending at 
the zero-coupling critical point remains. Thus only the zero-coupling point is left as a critical 
point in the parameter plane. 

Consider a self-similar sequence of parameters (A n ,c n ), at which the synchronous orbit 
of period 2 n has some given residues, in the U route for the global-coupling case. Rescaling 
the coupling parameter with the minimum scaling factor g ^ ^ (= 0.75), the sequence is 
transformed into a self-similar one for the N = 4 case of nearest-neighbor coupling. Hence, 
the critical behavior near the zero-coupling critical point becomes the same as that for the 
global-coupling case. 

The results for the nearest-neighbor coupling case with K = 1 extends to all the other 
nonglobal-coupling cases with 1 < K < y(^y^) for even (odd) N. For each nonglobal- 
coupling case with K > 1, we first consider a mode with index j min for which the scaling 
factor Sn ( K becomes the smallest one and the stability region U^(K, j min ) including a c = 
line segment. Here the value of j m { n varies depending on the range K. Like the K = 1 case, 
the zero c side of Un(K, j min ) including the c = line segemnt remains unchanged, whereas 
the other side becomes flattened by the bifurcation curves of the other modes with nonzero 
indices. Thus the overall shape of the stability diagram of the 2 n -periodic (n = 1, 2, 3, • • •) 
orbits born via synchronous supercritical PDB's becomes essentially the same as that for 
the nearest-neighbor coupling case. Consequently, only the U route ending at the zero- 
coupling critical point is left as a period-doubling route, and the critical behavior near the 
zero-coupling critical point is also the same as that for the global-coupling case. 



V. SUMMARY 

The critical behaviors of period doublings in the system of N symmetrically coupled 
inverted pendulums have been investigated by varying the two parameters A and c. As in 
the single inverted pendulum jy], the coupled inverted pendulums exhibit multiple period- 



doubling transitions to chaos with increasing A. We have studied the first three period- 
doubling transitions to chaos. For each period-doubling transition to chaos, it has been 
found that the critical behaviors vary depending on whether or not the coupling is global. 
For the global-coupling case the zero-coupling critical point and an infinity of critical line 
segments constitute the same critical set in the A — c plane, independently of N. However, 
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for any other nonglobal-coupling cases the structure of the critical set becomes different from 
that for the global-coupling case, because of a significant change in the stability diagram 
of 2 n -periodic orbits (n = 0, 1, 2, . . .). The critical scaling behaviors on the critical set have 
been also found to be the same as those for the abstract system of the coupled ID maps 

in. 
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TABLES 



TABLE I. For the case of the U route, the scaling factors [x n and v n in the scaling for the 
coupling parameter and the slope of the asynchronous residue at the zero-coupling critical point 
are shown in the second and third columns, respectively. 



n 




V n 


4 


-3.517 


-2.958 


5 


-2.904 


-2.627 


6 


-2.530 


-2.480 


7 


-2.495 


-2.522 



TABLE II. We followed, in the C route in Fig. 2(b), two self-similar sequences of parameters 
(A n , c n ), at which the pair of residues (i?o,n, Ri,n) of the synchronous orbit with period 2 n is (1, 0.1). 
They converge to both ends of the critical line segment. The scaling factors of the coupling paramter 
at the left and right ends are shown in the second and third columns, respectively. In both cases 
the scaling factors seem to converge to the same limit value fi ~ 2. 



n 






4 


3.66 


3.93 


5 


2.81 


3.04 


6 


2.02 


2.15 


7 


1.93 


1.99 



TABLE III. The scaling factors z^ n 's in the scaling for the slope of the asynchronous residue at 
the left and right ends of the critical line segment in Fig. 2(b) are shown in the second and third 
columns, respectively. 



n 


V n 




4 


2.528 


2.525 


5 


2.071 


2.072 


6 


2.001 


2.001 


7 


2.000 


2.000 
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FIGURES 



FIG. 1. (a) Stability diagram of the synchronous orbits of low period q = 1,2 in N lin- 
early-coupled inverted pendulums with the global coupling. Here A\ (= 0.575154 • • •) is just 
the first period-doubling transition point of the uncoupled inverted pendulum. The stable regions 
of the inverted stationary point, a symmetric 2-periodic orbit, and an asymmetric 2-periodic orbit 
are denoted by the SP, the SP2, and the ASP2, respectively. The horizontal (non-horizontal) solid 
and dashed boundary lines correspond to synchronous (asynchronous) PDB and PFB lines, respec- 
tively, (b) Phase portraits for A = 0.5. The phase flow of a symmetric 2-periodic orbit born via 
synchronous supercritical PDB is denoted by a solid curve, and its Poincare maps are represented 
by the solid circles, (c) Phase portraits for A = 0.57. The phase flows of a conjugate pair of 
asymmetric 2-periodic orbits born via synchronous supercritical PFB are shown: one is denoted 
by a solid curve, while the other one is denoted by a dashed curve. Their Poincare maps are also 
represented by the solid and open circles, respectively. 

FIG. 2. Stability diagram of synchronous asymmetric 2 n -periodic (n = 1,2,3,4) orbits of level 
n born via synchronous supercritical PDB's. PN denotes the stable region of an asymmetric orbit 
of period N (N= 2,4,8,16). The solid and dashed boundary lines represent the same as those in 
Fig. [I]. The stability diagram starting from the left (right) side of the ASP2 is shown in (a) [(b)]. 
Note its treelike structure. 

FIG. 3. Plots of the asynchronous residue R\^ n {A\,c) versus c near the zero-coupling critical 
point for n = 5, 6, 7. 

FIG. 4. (a) Plots of the asynchronous residue R\> n {A\,c) versus c near the critical line in 
Fig. 2(b) for n = 5, 6, 7. Their magnified views near the both ends q and c r are also given in (b) 
and (c), respectively. 

FIG. 5. Plot of the maximum asynchronous Lyapunov exponent u\ \ of the synchronous 
quasiperiodic orbit near the critical line in Fig. 2(b). This plot consists of 450 c values, each 
of which is obtained by iterating the Poincare map P 20 000 times to eliminate transients and then 
averaging over another 5000 iterations. The values of a% t i at both ends of the critical line are zero, 
which are denoted by solid circles. 

FIG. 6. (a) Second and (b) third stability diagrams of synchronous periodic orbits. Here A\ 
(= 3.829 784 • • •) and A\ (= 10.675 090) are just the second and third period-doubling transition 
points of the uncoupled inverted pendulum, respectively. The stable regions of the inverted sta- 
tionary point, an asymmetric orbit of period 1, a symmetric 2-periodic orbit, and an asymmetric 
N-periodic (N = 2,4,8) orbit are denoted by the IS, the ASP1, the SP2 and the PN, respectively. 
The solid and dashed boundary lines also represent the same as those in Fig. ||. 
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FIG. 7. First stability diagrams of synchronous periodic orbits near the c = line for the 
cases of (a) the quadratic and (b) cubic couplings. Here SP2 and PN (N= 2,4, 8) denote the stable 
regions of a symmetric orbit of period 2 and an asymmetric orbit with period N, respectively. 

FIG. 8. Stability diagram of synchronous periodic orbits in four linearly-coupled inverted 
pendulums with the nearest-neighbor coupling (K = 1). Each stable region is bounded by its solid 
boundary curves. For a synchronous orbit of period q, the PDB (PFB) curve of the mode with 

PD( PF} 

index j is denoted by a symbol g.- 
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